Topological quantum phase transition in bond-alternating spin- 1/2 Heisenberg chains 
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We investigate string correlations in an infinite-size spin-1/2 bond-alternating Heisenberg chain. By employ- 
ing the infinite matrix product state representation with the infinite time evolving block decimation method, 
a finite string correlation for extremely large lattice distances is directly observed, contrast to an extrapolated 
extreme value for finite size chains. We find that a topological quantum phase transition occurs between two 
different phases separated and characterized by two different long-range string orders in the space of bond- 
alternating interactions. Also, the critical exponent f} from the long-range string orders is obtained as /? = 1/12 
and the central charge at the critical point is obtained as c ^ 1 , which shows that the topological phase transition 
belongs to the Gussian universality class. In addition, it is shown that, for the topological quantum phase tran- 
sition, the phase boundary can be captured by the singular behavior of the von Neumann entropy and the pinch 
point of the fidelity per site. 
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I. INTRODUCTION 

Landau's symmetry breaking theory^ (LSBT) for under- 
standing the continuous phase transitions is one of the corner- 
stone theories in quantum many-body systems. Most phases 
and phase transitions could be understood by a local order pa- 
rameter characterizing a symmetry breaking Recently, a 
new type of phase transitions, for instance, appearing in the 
discovery of the quantum Hall effect"*, could not be under- 
stood well by local order parameters. Such quantum phase 
transitions seem to occur between two different phases with- 
out any explicit symmetry breaking. A non-local order, i.e., a 
so-called topological order^, rather than local orders, enables 
us to characterize the two phases. The quantum phase transi- 
tion is called topological quantum phase transitions (TQPTs) 
^il. Topological orders have been intensively studied in var- 
ious systems such as fractional quantum Hall state s'^''^ , Hal- 
dane phase chiral spin liquid s^i^ ^ Z2 spin liquids-*^--, and 
so on. 

In spin lattice systems, especially, a string order as a pro- 
totype example of non-local long-range order was introduced 
by Nijs and Rommelse^^, and Tasaki^"* and characterizes the 
Haldane phase. Since then, non-local string orders have 
been extensively applied to study quantum phase transitions 
in various spin systems such as anisotropic spin-1 Heisen- 
berg chains^^, frustrated antiferromagnetic Heisenberg spin 
chains^^, spin- 1/2^^'^** and spin-l^** alternating Heisenberg 
chains, XXZ spin-1 chain^, spin ladders^''^^, spin tubes"*^, 
and Kitaev spin mode l"''^ , and so on. Also, a TQPT has been 
shown between two phases respectively characterized by dif- 
ferent non-local string (topological) order parameters in Ki- 
taev spin models"''-. 

In this paper, we study a TQPT in an infinite spin-1/2 bond- 
alternating Heisenberg chain by introducing two non-local 
string orders. To compute the two long-range string orders, 
we employ the infinite matrix product state (iMPS) represen- 
tation-"*'^^ with the infinite time evolving block decimation 
(iTEBD) method developed by Vidal^^. Compared to finite- 
size lattice systems where a long-range order determined by 
an extrapolated value of correlations, the iMPS approach al- 
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FIG. 1: (Color online) Ground state phase diagram for the spin-1/2 
bond-alternating Heisenberg chain in the plane of the bond alterna- 
tions J and J' . The even- and odd-Haldane phases are denoted by 
yellow and green, respectively. The ferromagnetic phase occurs in 
the blue region. Note that the red solid line between the two Hal- 
dane phases denotes a continuous phase transition corresponding to 
the topological phase transition, while the red dashed lines indicate a 
discontinuous phase transition. 



lows us to calculate non-local long-range orders directly. 
The string correlations and the dimer correlations directly cal- 
culated from the iMPS groundstate show clearly that there are 
two topologically ordered phases characterized by two long- 
range string orders. From the topological characterization, we 
find that there are three phases including the ferromagnetic 
phase in the spin-1/2 bond-alternating Heisenberg chain (See 
Fig. [TJ. Furthermore, the central charge c ^ lat the phase 
boundary between the two topologically ordered phases and 
the critical exponent of the string orders = 1/12 show that 
the TQPT belongs to Gussian-type phase transition. We dis- 
cuss that the von Neumann entropy and the fidelity per lattice 
site (FLS) can capture the TQPT. 

This paper is organized as follows. In Sec. II, we introduce 
the spin-1/2 alternating Heisenberg chain model and discuss 
its groundstate energy. Section III devotes to discuss the topo- 
logical phase transition by means of the comparison between 
the odd-/even-string correlations and the dimer correlations. 



2 







0.8 0.9 1.0 1.1 1.2 
J7J 



O.odd 



O.even 





JVJ 



5> 



o 

LU 
I 



o 
LU 

CO 



4 ■ 



^ 0.1 

to 

o 

4^ 0.01 

in 0.001 



1e-04 



+ Data 
Fitting 



0.001 0.01 0.1 1 

5 



a Data 
— Fitting( (ln5-ln5o)/a ) 



0.01 



0.1 



FIG. 2: (Color online) Ground state energy per site eg and the odd- 
/even-bond energies eo,cv«i/orfrf as a function of J' with 7 = 1 . In the 
inset, the second order derivative of energy per site is plotted around 
the critical point J' = J = \. 



FIG. 3: (Color online) Estimation of the constants a and 6% from 
IMPS data. In the inset, Eq(Q) - Eo(6) as a function of 6 is shown. 



In Sec. IV, the phase transitions are discussed from the von 
Neumann entropy. The TQPT is classified based on the cen- 
tral charge via the finite-entanglement scaling. The FLS is 
discussed to detect the TQPT in Sec. V. Finally, our conclu- 
sion is given in Sec. VI. 



II. SPIN- 1 /2 ALTERNATING HEISENBERG CHAIN AND 
GROUNDSTATE ENERGY 

The spin- 1/2 alternating Heisenberg chains are described 
by the Hamiltonian^^--'^'-'^ 



H - ^^(JS 2i-l ■52/ + J'^ 2/ ■ ^ 2i+l), 



(1) 



where ^ - {S'\S^\S^) are the spin-1/2 operators. J and J' 
are two alternative nearest-neighbor exchange couplings. In 
order to cover the whole range of the parameter J and J', we 
take J - cos 9, J' = sin 6. For - n/A, this model is reduced 
to the conventional antiferromagnetic (AF) isotropic Heisen- 
berg chain with being gapless, while for 6 - or njl, the 
system will be characterized by decoupled singlets. When 9 
approaches the limit - -njl from - Q, Eq. ([TJ can be re- 
garded effectively as a 5 = 1 AF isotropic Heisenberg chain, 
which can be characterized by a gapful Haldane phase. 

From our IMPS groundstate in the Hamiltonian of Eq. ([T]|, 
one can calculate the groundstate energy. In Fig|2] the ground- 
state energy per site eo, as an average value of the energies for 
odd bond e() odd and even bond eo,eie«, is plotted as a func- 
tion of the alternation rate J' I J. We set the antiferromagnetic 
coupling J - \ and change J' . In the inset, the second-order 
derivative of the groundstate energy exhibits a singular behav- 
ior as J' I J approaches 1, which implies that the system under- 
goes a continuous (a second-order) phase transition across the 
critical point J' jj - 1 . 

To understand more the critical behavior of the alternat- 
ing Heisenberg chain in the view of how a bond-dimerized 



Heisenberg chain approaches the uniform limit {]' - J) in our 
IMPS approach, we also calculate the groundstate by using the 
Hamiltonian// = JY,i{\-S)S2i-v^2i+{'^+S)S2rS 2/+1, written 
in terms of the dimerization parameter 5. Here, J - (J + J')/2 
and 6^ (J'- J) /(J' + J) with < 5 < 1 . In the TableU a com- 
parison between the groundstate energies from DMRG^- 
and our IMPS approach is made for J - (J + J')/2 = 1 and 
the truncation dimension = 56. Note that our IMPS ground- 
state energies agree very well with those from DMRG. 

By virtues of a field theoretical approach. Black and 
Emery^'' have related the critical behavior to the four-state 
Potts model and found that the energy per site follows a 
power law with a logarithmic correction, i.e., Ei)(Q) - Eo(6) - 
a6'^^^/(ln6 - Indo) with the numerical fitting constants a ^ 
-2.2 and 6o ~ 110. By using DMRG, Papenbrock et alM 
have also found that a power law (oc 5^ '*^) is quite well in the 
range of 0.008 < 5 < 0.1, while logarithmic correction is need 
for small 5. From our IMPS representation, in Fig. [3] we plot 
the energy per site Eq{S) with the dimerization parameter 6 
for / = 1 . Eo{Q) - 1 /4 - In 2 is corresponding to the uniform 
limit. Figure [3] shows the power law behavior with the loga- 



TABLE I: Ground state energy per site Eq 



Dimerization 6 


DMRG £o 


iTEBD Eq(x = 56) 


0.000 


-0.443147'- 


-0.443146 


0.001 


-0.443166 


-0.443165 


0.002 


-0.443196 


-0.443197 


0.005 


-0.443333 


-0.443333 


0.010 


-0.443655 


-0.443655 


0.020 


-0.444537 


-0.444537 


0.050 


-0.448374 


-0.448376 


0.100 


-0.457246 


-0.457246 


1/3 


-0.517954 


-0.517954 


"Ref.m 






*Exact solution Eq = 


- In 2 + 1/4 given 


by Bethe Ansatz. 
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FIG. 4: (Color online) String order parameters Osir.odd/even and ferro- 
magnetic order parameter O^m as a function of 9 with J = cos 9 and 
y = sine. 



rithmic correction, £0(0) - £'o(^) - a6*^^/(ln6 - \n6o) with 
the numerical fitting constants a w -2.21 and 60 ~ 1 14. 

III. TOPOLOGICAL QUANTUM PHASE TRANSITION 

In the antiferromgantic Heisenberg Hamiltonian of Eq. ([TJ 
with the bond alternation, the transition between the two 
dimerized states at J' /J = 1 does not involve any change of 
symmetry although the second derivative of the groundstate 
energy exhibits its singluar behavior indicating phase transi- 
tion in Fig. |3] The absence of any explicit change of sym- 
metry implies no local order parameters characterizing each 
phase. However, it manifests as a change of topological order, 
i.e., string order To see this, let us introduce two string order 
parameters based on the bond alternation: 



2j-2 



^ srr,odd 



lim 

l'-;l->" 



-4<5^,exp[/;r J] 

k=2i+l 
2j-\ 

-4<5^,,iexp[/7r S'^jS^.p 



k=2i+2 



(2a) 



(2b) 



where a - x,y, and z- Due to the SU(2) rotationally invariant 
ground states, (9j,^ = (9*,,. = (9~j„.. From our IMPS ground- 
wavefunction, we directly calculate the defined string orders. 
In Fig. m the order parameters are displayed as a function of 
the angle variable 6 = tan"' J' /J. As we mentioned, the string 
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FIG. 5: (Color online) String order parameter Os,r,odii as a function 
oiJ'IJ with J = 1. In the inset, the string order parameter O^tr.odd is 
plotted as a function of 5''*. The straight line, as a guide of eye, is 
also drawn. 



orders are the saturated values of string correlations for a large 
lattice distance. Figure|4]clearly shows that the odd string or- 
der is finite in the -7t/2 < 9 < n/A phase and zero otherwise, 
while the even string order is finite in the ;7r/4 < < n. These 
results imply that each system parameter range can be char- 
acterized by each long-range string order Also, as the system 
parameter crosses the critical point J' = J, the system under- 
goes a TQPT between the two topologically ordered phases. 
Similar to the spin-1 Heisenberg chain understood by the hid- 
den Z2 ® Z2 breaking symmetry^-' '*'^, a similar hidden symme- 
try breaking may occur for each phase. Further, we plot the 
string order parameter Ostr.odd a function of J' with J - 1 
in Fig. |5] In the inset of Fig. |5] the string order parameter 
Osir.odd is plotted as a function of the dimerization parameter 
6. Note that the string order parameter Ostr.odd scales as 6^^^ 
— . As a result, the critical exponent is given as /? = 1/12 via 
Ostr.odd °^ Hence, the TQPT is in the same universality 

class with the Gaussian phase transition. 

When the 9 approaches -;7r/2 or n, the system can be re- 
garded as an effective spin-1 antiferromagentic Heisenberg 
chain because of the strong odd- or even-bond dimeriza- 
tions. In these limits, the numerical values of our string 
orders Ostr.odd/even = 0.3873 become very close to the re- 
sults, for the antiferromagentic spin-1 Heisenberg chain, 
(9f,=' = 0.37434447 from the IMPS approach^ and 0?=.' = 
0.37432509 from the DMRG method"*'. For -n < < -n/l 
(J' < and J < 0), both the odd and even string orders be- 
come zero, while the ferromagnetic local order OfM = {Si) 
become finite, which indicates that the system is in the fer- 
romagnetic phase. In the magnetization in Fig. [H one may 
notice a non-zero value of the local magnetization (Si) at 
- 7r/4, i.e., an isotropic Heisenberg chain that should have 
zero magnetization. The non-zero value of the magnetization 
at - 7r/4 is due to finite truncation dimension However, 
it scales down to zero in the limit of infinite truncation dimen- 
sion, i.e., the thermodynamic limit—. 

Actually, Hida^^ introduced a similar string order param- 
eter to argue a crossover between a Haldane-gap phase and 
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FIG. 6: (Color online) Dimer orders Oo.even as a function of 9 with 
J = cos 6 and J' = sin 0. 

a 'dimer phase' in the spin- 1/2 bond alternating Heisenberg 
chain for J > J'. However, his 'dimer phase' is originated 
from the symmetry of the Hamiltonian rather than a broken 
symmetry of groundstate because both the Hamiltonian and its 
groundstate possess the same symmetry, i.e., two-site transla- 
tional invariance. In this sense, Hida has found the existence 
of a topological order for J' > J in the spin- 1/2 bond alter- 
nating Heisenberg chain for J > J'. However, for J > J', 
no crossover happens in the system because his 'dimer phase' 
might be rather a 'dimerized state' due to the bond alterna- 
tion. To make this point clearer, we calculate the odd- and 
even- dimer orders defined by 

lim <Oo(20Oo(2;)>, (3a) 

OD,odd^ lim <(9o(2/ + 1)0d(2;)>, (3b) 

where ^^(O - 5*,- ■ ^i+i - 5*,+i ■ In our system, we cal- 
culate the even- and odd- dimer orders Oo.odd - -OD,even- In 
Fig. |6] the even-dimer order is displayed as a function of the 
angle variable 0. Both the odd- and even-dimer orders are fi- 
nite in the two Haldane phases, which implies that the two 
Haldane phases cannot be distinguished by the dimer orders 
although both the odd- and even-dimer orders disappear at the 
critical point - njA. These disappearances of the dimer 
orders at the critical point are originated from the symmetry 
of the Hamltionian because the Hamiltonian is one-site trans- 
lational invariance and its groundstate has the same symme- 
try of the one-site translation invariance. Consequently, the 
dimer orders are not the order parameters characterizing the 
two Haldane phases although the dimer orders exist in the two 
Haldane phases. Hence, the two long-range string orders char- 
acterize the two Haldane phases, respectively, and the TQPT 
occurs between the two Haldane phases. 

IV. VON NEUMANN ENTROPY AND CENTRAL CHARGE 

At a quantum critical point, characteristic common singu- 
lar behaviors of thermodynamics properties, depending only 
on few features such as dimensionality and symmetry, allow 
us to classify quantum phase transitions by using the concept 
of the universality classes, i.e., a type of quantum phase tran- 
sition. Especially, the central charge'*^ can be used to clas- 
sify a universality class of the TQPT at J' /J = 1. In the 



iMPS representation, at a quantum critical point, one can ob- 
tain the central charge c via a so-called finite-entanglement 
scaling exponent k numerically from the unique behaviors of 
the correlation length ^ = ax" and the von Neumann entropy 
S = (c/f/6) logj;^^'*. Here, a and x are a numerical fitting 
constant and the truncation dimension, respectively. Thus, in 
order to get more insight on the TQPT, we will discuss the von 
Neumann entropy and the central charge ai J' jj - 1 in this 
section. 

Actually, the von Neumann entropy can quantify quantum 
entanglement as a good measure of bipartite entanglement be- 
tween two subsystems of a pure state"*^-"*^. In our iMPS repre- 
sentation, the spin chain can be partitioned into the two parts 
denoted by the left semi-infinite chain L and the right semi- 
infinite chain R. The von Neumann entropy is defined as 
S - -Irgi log2 Ql - -Trg« log2 qr in terms of the reduced 
density matrix qi or pR of the subsystems L and R. The von 
Neumann entropy for the semi-infinite chains L or R can be 
expressed as 

Si--Y,^U°S2 4a^ (4) 

where /l,;u,'s are diagonal elements of the matrix A that could 
be directly obtained in the iMPS representation. This is 
because, when one partitions the two semi-infinite chains 
L(-oo, ■ ■ ■ , /) and R(i + 1, ■ ■ ■ , oo), one gets the Schmidt de- 
composition I^F) = 2^=1 ^i,ffl0L)l0fl)- From the spectral de- 
composition, are actually eigenvalues of the reduced den- 
sity matrices for the two semi-infinite chains L and R. In our 
two-site translational invariant iMPS representation, there are 
two Schmidt coefficient matrices Aa and Ag that describe two 
possible ways of the partitions, i.e., one is on the odd sites, 
the other is on the even sites. From the Aa and Ag, one can 
obtain the two von Neumann entropies depending on the odd- 
or even-site partitions. 

In Fig.|2](a), we plot the von Neumann entropies Sodd and 
Seven as a function of the angle variable 6 = tan"' J' /J. Both 
the entropies for odd and even bonds show their singular be- 
haviors at the same values = -n/l, jijA, and ji. Note that the 
singular behaviors of the von Neumann entropies correspond 
to the quantum phase transition points from the order param- 
eters in Fig. |4] Hence, the von Neumann entropies give the 
same phase diagram from the order parameters. 

In Fig. |2](b) and (c), the correlation length ^ and the von 
Neumann entropy S are plotted as a function of the trunca- 
tion dimension at the critical point 9 - tt/A {J' I J - 1). 
The truncation dimensions x are taken from 16 to 64. It is 
shown that both the correlation length ^ and the von Neumann 
entropy S diverge as the truncation dimension x increases. 
From a power-law fitting on the correlation length ^, we have 
K - 1.336 and a - 0.255. As shown in Fig.|7](c), our numeri- 
cal result demonstrates a linear scaling behavior, which gives 
a central charge c ^ 1.004 with k - 1.336. Our central charge 
is close to the exact value c - \ predicted by conform field 
theory for an isotropic antiferromagnetic Heisenberg chain. 
Consequently, the TQPT between the odd- and the even- Hal- 
dane phases at J' /J - 1 is a Gaussian transition which is 
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FIG. 7: (Color online) (a) Von Neumann entropies S „f„ wAS odd as a 
function of the parameter 6. The von Neumann entropies for odd- and 
even-bonds show a peak at 6 = 7r/4 and discontinue at 6 = -nil, n. 
(b) Correlation length f as a function of the truncation dimension^ at 
the critical point 9 = n/4. The power curve fitting £ = ax^ yields a = 
0.255 and k = 1.336. (c) Scaling of the von Neumann entropy S with 
the truncation dimension at the critical point 9 = n/4, i.e., J' = J 
corresponding to the conventional spin- 1/2 AF isotropic Heisenberg 
chain. For k from (b), the linear fitting S = (ck/6) log2X + b yields 
b = 0.369 and the central charge c a 1.004. 

characterized by a central charge c - 1 and the occurrence 
of a phase transition between two gapful phases. This result is 
consistent with the classification of the TQPT from the critical 
exponent /? = 1/12. 

V. FIDELITY PER LATTICE SITE 

As we discussed in the previous section, the von Neumann 
entropy, as an entanglement measure, can detect the topolog- 
ical phase transition at J' /J - 1. In this section, we study 
a fidelity per site as another universal indicator for quantum 
phase transition. From our iMPS groundstate wavefunction 
1^(6)} in terms of the angle variable 6 - tan"' J' /J, one can 
define a fidehty F(6', 6) = \{'V{ff)\'V{6))\ between the ground- 
state wavefunctions for two different control parameters 9 and 
ff . A fidelity per lattice site (FLS)^ d can be defined as 

\nF{e',e) 

\nd{e',e) = lim — -, (5) 

where N is the system size. In Fig. |8] the FLS d{ff , ff) is plot- 
ted in the two dimensional control parameter space. The three 




FIG. 8; (Color online) Fidelity per site d(9',6) surface as a function 
of two parameter 9' and 9 for a bond alternating Heisenberg chain 
with J = cos 9 and J' = sin 9. 



singular points (ffyO) - (7r/4,7r/4), {-nil, -nil) and {7i,n) 
are observed on the FLS surface, which indicates that there oc- 
cur quantum phase transitions when control parameters across 
these values. The continuous behavior of the FLS function 
across the "pinch" point {ff ,0) - {7tlA,nlA) implies that a 
continuous quantum phase transition occurs for the topologi- 
cal phase transition. The discontinuous behaviors of the FLS 
function at two other points are corresponding to the first- 
order quantum phase transitions"*^. The pinch points of the 
FLS correspond to the quantum phase transition points from 
the order parameters in Fig. |4]and the von Neumann entropies 
in Fig. Ha). The FLS gives the same phase diagram from the 
order parameters and the von Neumann entropy. Hence, it is 
shown that a TQPT can be detected by FLS. 



VI. CONCLUSION 

Quantum phase transitions have been investigated system- 
atically in a spin- 1/2 bond-alternating Heisenberg chain by 
using the iMPS representation with the iTEBD method. By 
calculating the odd- and even-string orders and the feiTomag- 
netic order, the three phases, i.e., the odd- and even-Haldane 
phases and the feiTomagnetic phases, were found in the plane 
of the bond-alternating interactions. The TQPT between the 
odd- and even-Haldane phases was classified as a Guassian- 
type phase transition from the central charge at the critical 
point and the critical exponent of the string orders. Also, it 
was clearly shown that the FLS and the von Neumann entropy 
can detect the TQPT. 
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